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Abstract 



We derive the semiclassical series for the partition function of a one- 
dimensional quantum-mechanical system consisting of a particle in a single- 
well potential. We do this by applying the method of steepest descent to 
the path-integral representation of the partition function, and we present a 
systematic procedure to generate the terms of the series using the minima 
of the Euclidean action as the only input. For the particular case of a quar- 
tic anharmonic oscillator, we compute the first two terms of the series, and 
investigate their high and low temperature limits. We also exhibit the non- 
perturbative character of the terms, as each corresponds to sums over infinite 
subsets of perturbative graphs. We illustrate the power of such resummations 
by extracting from the first term an accurate nonperturbative estimate of 
the ground-state energy of the system and a curve for the specific heat. We 
conclude by pointing out possible extensions of our results which include field 
theories with spherically symmetric classical solutions. 
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I. INTRODUCTION 



Semiclassical series were introduced in Quantum Mechanics by the pioneering works of 
Brillouin [|IJ , Kramers |2| and Wentzel || . In order to solve the time-independent Schrodinger 
equation for slowly- varying potentials, they proposed an ansatz for the wavefunction of the 
form (here particularized to one dimension) 

^( x ) = e iS (*)A (1) 

that led to a Ricatti equation for S(x). The assumption of slow variation of the potential 
allowed for an iterative procedure in solving that equation. As a result, one obtained a series 
expansion of S(x) in powers of h, i.e., S(x) = J2™=o S n (x)h n , with S n+ i(x) given in terms of 
S n (x) by means of a recursion relation, and S (x) satisfying the Hamilton- Jacobi equation 
for a particle with the same energy and external potential as those of the Schrodinger 
equation. Thus, the zeroth order approximation could be associated to a classical path, 
from an arbitrary x to x, whose classical action Sq(x) = S c i(xo — > x), apart from a sign, 
should be used in ([I]). The choice of xq just fixed a normalization constant. 

The terms of the series derived in references |]-[| were plagued with divergences that 
occured at the turning points of the classical motion, a problem that came to be known as 
the connection problem. Dunham bypassed this difficulty by turning x into a complex 
variable z, and using different linear combinations of e ±tS<yX ^ h as asymptotic series for ij)(x) 
along the real axis. Assuming a potential with only two classical turning points x a and Xb 
(x a < Xb), three different linear combinations should be used as x < x a , x a < x < Xb or 
x > Xb- Demanding that ip(z) be real, bounded and single- valued on the real axis, and 
writing S(z) = J z p(z') dz', he could arrive at the quantization condition 

p(z)dz = 2im% (n = l,2, . . .), (2) 

for any closed contour encircling the segment (x a , Xb)- The semiclassical series for ip(z) led to 
semiclassical series for S(z) and p(z) which, together with (g), provided an expression for the 
energy levels as a power series in h, generalizing the Bohr-Sommerfeld condition. That series 
was used as a starting point for the works of Bender, Olaussen and Wang Q and Balian, 
Parisi and Voros 0,0], who refined such asymptotic representations to obtain highly accurate 
estimates of the energy levels of the single-well quartic anharmonic oscillator. These could 



be compared to the estimates obtained by Bender and Wu PHT0| from the investigation of 
the large order behavior of perturbation theory jTT-14 . 



Apart from energy eigenf unctions and eigenvalues, semiclassical series were also written 
down for Green's functions by Balian and Bloch [fT5 |. Defining the resolvent operator G 
through 

(H-z)G = 1 (3) 

for a complex energy z, H being the Hamiltonian of the system, they used an ansatz of the 
form 

G{x, x'; z) = A{x, x') e s{x > x '' z)/h . (4) 



2 



From the defining equation (^j, under assumptions similar to those in the first paragraph, 
they started from zeroth order estimates S (x, x'] z) and A (x,x') by neglecting nonleading 
terms in h. They found that Sq satisfied a Hamilton- Jacobi equation for complex energy 
z, solved by a classical trajectory from x to x' . The amplitude A (x,x') turned out to 
be related to derivatives of the classical So- Just as before, the zeroth estimates could 
be used to generate a whole series expansion by an iterative procedure. Note that z = 
E + £7, the analytic continuation to complex energies being used to circumvent problems 
with singularities (turning points in one dimension; caustics in higher dimensions). 

In all the discussions mentioned, the first term of the semiclassical series was obtained 
from quantities related to some classical trajectory. Path integral representations for cor- 
relation functions 



16-20 



which sum over trajectories, should therefore provide a natural 
setting to derive semiclassical series. Indeed, approximating quantum-mechanical path inte- 
grals by the stationary phase method leads to zeroth order estimates given by the paths that 
solve the classical equations of motion, the saddle points of the action functional. Thus, 



(z|e-4<*-*V) 



x(t)=x 
x{t')=x' 



[Dx{t)\ e iS[x{t)]/h , 



(5) 



may be time Fourier transformed to yield The saddle-point contribution, and fluctua- 
tions around it, will yield the zeroth order estimates Sq and A$ mentioned before. 

Although semiclassical quantization by means of path integrals became widely used, as 
evidenced by Gutzwiller's work and by extensions to field theories by Dashen, Hass- 



lacher and Neveu [p^ -pql, almost all discussions never went beyond the first term of a semi- 
classical series. Notable exceptions were the works of DeWitt-Morette PB" |, for arbitrary 



potentials in Quantum Mechanics, and Mizrahi |p7[| , for the single- well quartic anharmonic 
oscillator in Quantum Mechanics (more recently , a derivation of the higher order terms 
of the series using transport-type recurrence equations also became available). However, 
these contributions to the mathematical physics literature did not receive the attention they 
certainly deserved in more applied work. We were no exception, and only became aware of 
these articles after we had rederived the semiclassical series using the methods of section ||. 



Semiclassical methods for finite temperature field theories [E9H31j] also remained re- 



stricted to derivations of the first term of a semiclassical series [32], even when the problem 



was reduced to Quantum Statistical Mechanics |3"3| , |3"3|] , viewed as field theory at a point (zero 
spatial dimension). Some references resorted to extensions to the complex plane |35|-p7| to 
include complex paths required to describe Fourier transformed quantities but, again, those 
treatments were not concerned with obtaining the whole series. 

The central objective of the present article is to bridge this gap. Thus, we undertook 
the task of finding a systematic path integral procedure to generate a semiclassical series for 
Quantum Statistical Mechanics which led us to the construction of each term of the series 
from the knowledge of the solution(s) of the classical equations of motion. We concentrated 
our attention on the partition function, for which use of the method of steepest descent 

m 



only required real solutions as saddle-points 



The restriction of our analysis to one- 



dimensional quantum- mechanical systems (i.e., scalar field theories at a point and at finite 
temperature) allowed us to construct the semiclassical propagator needed to generate the 
terms of the series. 
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This article is organized as follows: section || presents the derivation of the semiclassical 
series for a generic potential of the single-well type, both in quantum-mechanical language 
and in field-theoretic language, the latter allowing for a simple connection with the works 
of references p6|j27| . It should be remarked, however, that our presentation is quite simple, 
being a natural extension of textbook material ||14j|, and having profitted greatly from the 



clear account of reference [ 37j • To emphasize that our construction is free of turning-point 
singularities, details of the derivation were worked out in Appendices |A] and [BJ section |TJ 
discusses the single-well quartic anharmonic oscillator. There, we constructed the two first 
terms of the series explicitly, and looked at relevant limits. We also computed the ground- 
state energy and the specific heat, as illustrations. Appendices y, |B] and [E] complement the 
calculations in the text; section [TV] presents our conclusions, points out directions for future 
work, and lists a number of situations where our results could be applied. 



II. STATISTICAL MECHANICS 



A. Quantum- mechanical path integrals 



The partition function for a one-dimensional quantum-mechanical system consisting of a 
particle of mass m in the presence of a potential V(x) in equilibrium at inverse temperature 
/3 can be written as a path integral: 

poo px(f3h)=x 

Z((3) = dxj [Vx{r)\ e- s ' r \ (6) 

J~oo Jx(0)=xq 



S[x] 



dr 



Hi +v[x) 



(7) 



For convenience we define the dimensionless quantities q = x/x^, 9 = u N r, 9 = Phu^, 
U(q) = V(x]yq) /moj 2 N x 2 N and g = Ti/muj N x 2 N , where u^ 1 and xn are natural time and length 



scales of the problem, respectively, 
function as 



Z(G) 



In terms of these quantities we rewrite the partition 



r q(0)=q o 

/o / , , Pq{0)] e 



1/9 



(8) 



(9) 

where the dot denotes differentiation in 9. 

We generate a semiclassical series for Z(<d) by: (i) finding the minima q c {9) of the 
Euclidean action J, i.e., the stable classical paths that solve the Euler-Lagrange equation 
of motion, subject to the boundary conditions; (ii) expanding the Euclidean action around 
these classical paths; (iii) deriving a quadratic semiclassical propagator by neglecting terms 
higher than second order in the expansion; (iv) using that propagator to compute higher 
(than quadratic) order contributions perturbatively. 



d9 



^ 2 + U(q) 
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For the sake of simplicity, we shall restrict our analysis to potentials of the single-well 
type, twice differentiable, and such that U'(q) = only at the minimum of U, which we shall 
assume to be at the origin (see Fig. [I]). This guarantees that, given q and 0, there will be 
a unique classical path satisfying the boundary conditions. Multiple-well potentials force us 
to consider more than one classical path for certain choices of go and 0. This phenomenon 
has been analyzed, for a double-well type potential, using the language of catastrophes and 
bifurcations p9||. Semiclassical series for the double- well quartic oscillator will be presented 



elsewhere 38 



The Euler-Lagrange equation (U' = dll/dq) 

q-U'(q) = 0, (10) 

subject to the boundary conditions q(0) = q(Q) = qo, describes the motion of a particle in 
the potential minus U . Its first integral is 

1 -q 2 = U(q)-U(q t ) } (11) 

where q t denotes the single turning point (since we have an inverted single well) of the 
motion, defined implicitly by 

9-2 r -At. (12) 



'/() 



v(q,q t ) 



where v(q, q') = sign(g' — q)y 2[U(q) — U(q')\, and equation flT2| ) is a consequence of integrat- 
ing (pT|). Thus, for a single well, given q and 0, the classical path will go from q , at 9 = 0, 
to q t = qt(qo, 0), at 9 = 0/2, and return to q at 9 = 0. (Note that sign(g 4 ) = sign(g )-) 
The action for this classical path has a simple expression in terms of its turning point: 



I[q c ] = eU(q t ) + 2 / dqv(q,q t ), (13) 



where we have used flTT|). The first term in ( P^|) corresponds to the high-temperature limit of 
Z(Q), where classical paths collapse to a point (q t — * g ). The last term will be negligible for 



potentials that vary little over a thermal wavelength A = h^J p/m. However, by decreasing 
the temperature it will become important and bring in quantum effects. 

We now expand the action around the classical path. Letting q{9) = q c {9) + rj(9), with 
77(0) = 77(0) = 0, we obtain 



where 



I[q] = I[q c }+I 2 [ri} + tm, (14) 



UV] = \ [ d9 {rf(9) + U"[q c {9)\ V 2 (0)} , (15) 



r e 00 1 r e 

SI[ V }= d9 6U(9, V ) = Y,- d9U^[q c (9)}v n (0). (16) 

«/ q TXi > J 



n=3 U - J0 
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Inserting (|I4"D into (H) and expanding e in a power series yields 



Z(Q) 



r r)(&)=0 00 1 

dq oe - I[qc]/9 I Pv(0)} e- 1 ^ 9 £ 



r,(0)=0 



m=0 



(17) 



If we interchange, a la Feynman ||16|| , the integral over deviations from the classical path 
with the integrals over various times {0i}, coming from insertions of (|T6|) , we break up the 
path-integral into pieces which go from zero to the various {9i} and, finally, to 6. Thus, 



z(e) 



dq e 



-i[ic]/g 



00 (_l)n 

G c (o,o ; e,o) + 5:^^ r G n (go,e) 



(18) 



dr ]j ]G c {0,0;9 1 ,r h )5U{9 1 ,r h ) 



n / m j 

j= i J o 

xG c (6»i, 771; 2 , ^2) • • ■ <H7 (0 n , r? n ) G c (6 n , r] n ; 9, 0) 



Ji<0 2 <...<0„ 



(19) 



with G c (9i, rji, 9 2 , ?7 2 ), the semiclassical propagator, emerging naturally from doing the path- 
integral in each piece: 



G c {9i, T]x, 6 2 , T] 2 ) 



r\(Qi)=r\i _ 
Vifil)=Vl 



'[V V (9)]e 



-h[0i,02;v]/9 



(20) 



h[e^e 2 -r l ] = l -f d i d9 {rj 2 + U'% 



v 2 }. 



(21) 



It remains to show how one can obtain that propagator from the knowledge of the classical 
path. This is the central point in the whole procedure. For this, we use the fact that the 
action I 2 is quadratic in 77, and so the path integral in fl2"U|) is completely determined by the 
extremum rj e (9) of I 2 [9i, 9 2 ;r]}, which satisfies 



V-U"[q e (9)]V = 0, 
subject to the boundary conditions 11(81) = r/i and rj(9 2 ) = r\ 2 . Thus, 

G c (0i,77i;0 2 ,7fc) = G c (9 1 ,0;9 2 ,0)e- I ^ e2 ^, 
where, after an integration by parts, 



(22) 



(23) 



h[9i,e 2 ; Ve ] = -[r l2 r J e(82)-ViV e (0i)}- 



(24) 



We can obtain r] e (9) by finding the linear combination of any two linearly independent 
solutions, rj a {9) and rjb(9), of (|22|) which satisfies rj e (9i) = r\\ and rj e (9 2 ) = r\ 2 . The result is 



Ve(0) 



Vl n(9,9 2 )+r ]2 n(9 1 , l 

n(9 u 9 2 ) 



(25) 
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where 



We may then write 



[W 12 7]l + W 21 T)\ - (W U + W 22 ) T/ir/a], 



2tt 



12 



(26) 



(27) 



where fly = fl(9i,9j) and = dVtij/d9j. (Note that Wu is the Wronskian of r\ a and i] b 
computed at 6>j.) 

Explicit expressions for rj a (9) and rjb(9) can be obtained as follows. By differentiating 
(|10|) with respect to 9, one can verify that rj a (9) = q c {9) satisfies (p2|). For the second 
solution, we take rj b {9) = q c {9) Q(9), where Q{9) is defined as 

r 9 d9' 

Q(9)=Q(0)+ — — (28) 



q 



for 9 < 0/2, Q{9) = -Q(B - 9) for 9 > 0/2, and Q(0) is chosen so as to make r} b (6) 
continuous at 9 = 0/2 (see Appendix A). One can easily check, using (0), that rj b {9) indeed 
satisfies (|22"|). (Alternatively, one could use a procedure introduced by Cauchy |26],|27|], and 
differentiate the classical solution q c (9) with respect to any two parameters related to its 
two constants of integration.) We can now write explicit expressions for Q\ 2 and Wij\ 



n 12 = qc(0i)UO 2 ) [Q(9 2 )-Q(9 1 )), 



q c 



W tj = q c (°i) UWj)] [Q(0j) ~ Q(6i)] + - 



(29) 



(30) 



As a final step, the pre-factor in fl2"3"D is derived in Appendix ^, using the methods of 

1/2 



Refs. P^]57| . The result is 



G c (9i,0;9 2 ,0) 



2ngQ 



12 



(31) 



From (|30D , one easily finds Wu = 1. Therefore, our quadratic semiclassical propagator is 
given by 



Gc(9i,r]i,0 2 ,r] 2 ) 



1 



exp 



1 



2gtt 



(W 12 rjl + W 2l ril - 2 77x773) 



12 



As promised, it is completely determined by the classical solution. 

Finally, we note that the van Vleck determinant A is a by-product of 

A(g , ©) = G- 2 (0, 0; 0, 0) = 2ng fi(0, 0) = ing g c 2 (0) Q(0). 

Using (Tn|) and ( |A7[ ) one can express A as 

= Airg [U(q t ) - U(q )} ( dQ\ 

u f (qt) \dqt) qo - 



(32) 



(33) 



(34) 



Together with (13), this shows that one does not need to know q c (9) in order to write the 
first term in the semiclassical series (|TB); it is enough to know qt(qo, 0)- 



B. Field-theoretic approach 



We may generate the semiclassical series in an alternative way, which connects it quite 
naturally to the diagrammatics of field theory. 

The summation in (|T7|) can be written more explicitly as 



y J_ ( J_M ) " = i + v TT 



OO 1 



n,-=3 n J ! J ° 



(35) 



As a consequence, one is led to compute integrals of the following type: 

r r)(S)=0 
Jr;(0)=0 



(36) 



Such integrals emerge naturally as functional derivatives of the following generating func- 
tional: 



rr,(0)=O 

Z[J)= [D V (9)]e 

Jr?(0)=0 



-±{l2[0,&;r,}-f®dej(6) V {0)} 



(37) 



Indeed, 



5J(9 1 )---5J(9 k 



(38) 



j=o 



In order to compute Z[J), we define 



V(9)=V(9)+ / d9'g(9,9')J(9') 
Jo 



where r/(0) = 77(0) = 0, and 8(9,9') satisfies 

r 9 2 ^ 



d9 2 



+ tf"[g c (0)] 0(0, 0') = 5(9 - 9'), g(0, 9') = 0(6, 0') = 0. 



(39) 



(40) 



Inserting ( j3~9"D in fl37D , and noting that [Vr](9)] = \Vfj(9)), we obtain 

Jrj(0)=0 

The path integral in (f4~I|) has the same form as the one in (0), so we finally arrive at 



(41) 



Z[J] = G c (0,0;6,0) exp 



— [° d9 [ d9'J(9) Q(9, 9') J(9') 
2g Jo Jo 



(42) 



Using this result, we can now calculate fl38[). The result is simply 

• ■ ■ V(0k)} = 9 h/2 G c (0, 0; 9, 0) £ <?(0 4l , i2 ) ■ ■ • 0(0*^ A). ( 4 3) 



if k is even, and zero otherwise. J2p denotes sum over all possible pairings of the 9i r 
Inserting this into (|T7|) and (^) yields the semi classical series for Z(Q). 

We still have to solve Eq. (f40|). This can be easily done if one notes that, for 9 ^ 8', it 
has the same form as Eq. fl22|) . Therefore, Q(9,9') can be constructed, just as rj e {9) itself, 
as a linear combination of the rj a {9) and rjb(9) defined previously: 



a-r ]a (0)+b-7i h (0), 9<9' 
a+rj^+b+r]^), 6>6'. 



(44) 



Continuity imposes 



g(6' + e,6') = g(6'-e,6'), 



(45) 



whereas (HDD leads to 



d 

89 v ' ' 



(46) 



with e — > + . (|4q) , ( [46|) and the boundary conditions completely determine the coefficients 
in (H5). The final result is 



n(O,0<)n(0>,e) 

n(o,e) ' 



(47) 



where #<(#>) = min(max){6 1 , 6 1 '}, and fi(#i,#2) is the function defined in 



III. THE SINGLE- WELL QUARTIC OSCILLATOR 

In this section, we study the potential 

V(x) = - muo 2 x 2 + - Ax 4 . 



(48) 



Choosing = ui and xjy = \Jmuj 2 /\, and introducing the dimensionless quantities of 
section H, we have g = Xh/m 2 uj 3 and 



Integrating ( |llf) leads to W\ W^ 



q c {9) = q t nc(u e ,k), 



(49) 



(50) 



where nc(u, k) = l/cn(w, fc) is one of the Jacobian Elliptic functions [40-42 , and 

k 



u e = yl + qt 



0' 
2J> 



2 + q? 



\Hl + q?Y 



(51) 
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For future use, we note that (|51| ) can be rewritten as 

26 -6 , , /2(1 - k 2 ) 



Uq 



2 V / 2F^I' m V 2P-1 ' 
The relation between go and qt is obtained by taking 9 = in (]50|) : 

?o = <?c(@) = ft ncM Q . 

(We shall often omit the fc-dependence in the Jacobian Elliptic functions.) 
The action for the classical path ( |50|) is 

I[Qc] = 6 U(q t ) + V2 [ M dq v/(g 2 + g f 2 + 2)(g2-g|). 

J\qt\ 

Performing the integral [Ref. |5|, formula 3.155.6] and replacing go by the r.h.s. of 
obtain 



(52) 



(53) 



(54) 



we 



%c] = e(]g t 2 + ig t 4 )+^-Vl + S? 



+ sntte (l + ^g 2 nc 2 w e ) yi + t;<? 2 (1 + nc 2 M 6 )| , 



(55) 



where E(<£>, fc) denotes the Elliptic Integral of the Second Kind and tpg = arccos[g c (#)/go] 
arccos(cnMe). 

For the construction of the quadratic semiclassical propagator we shall need 



T} a {0) = q c (0) = q t V 1 + qfsnue dnu nc u e 



(56) 



and 



g(0)=g,- 2 (l + g 2 )- 3 / 2 

cn Ug dn Ug 



+ (k 2 ~ 1) 



cn ug sn Ug 



sn -u e dn ug 

We may then obtain rj h {6) = q c {6) Q{6) and, thus, Q i2 and W12 from (|29p and (|3 



use of (|32|) will yield the desired propagator . 

For the series expansion of the partition function ([18]), we shall need 



5U(9,r ] ) = q c (6)r ] 3 + -r ] 4 , 



(57) 
Finally, 

(58) 



obtained from ([IB]). Therefore, we have to consider not only the usual quartic vertex, 
but an additional time(#)-dependent cubic term. This completes the set of ingredients 
needed to write down a semiclassical series for any correlation. In the next subsection, we 



shall concentrate on the first term of the series (0) for Z(Q), which yields the quadratic 
approximation. 
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A. The quadratic approximation for Z(Q) 



From the knowledge of the classical action and of the Van Vleck determinant, we define 

/oo 
dq e' IM/g A~ 1/2 (59) 
-oo 

as the quadratic approximation to Z(Q). To perform the integral over go one must write 
I[q c ] and A solely in terms of go (and G), but except in rare cases this is not an easy task. 
Usually, it is much simpler to write these quantities in terms of q t [see Eq. (|53D and Appendix 
|Cj, and so it is natural to trade go for q t as the integration variable in (|59|) . This is much 
simplified by the fact that the Jacobian of the map go — > qt is simply related to the van 
Vleck determinant. In fact, Eqs. (|T2|) and ([34]) imply 

/«9g \ _ (dQ/dqtU 1 , Jd®\ _ U'(q t )A 

[WJe ~ ~(dQ/dq ) qt ~ 2 V ^ qt) [dq- t ) qo ~ ~^ 9 v(q ,q t y (60) 

Eq. ( |59"D then becomes 

i r« ^ u\ qt ) a 1 / 2 __ I[ftJ/ , ,«$ ^ 



3,(6) = -— / dq t y ; l > e- 1 ^ 9 = / dq t D{q t ,Q)e- 1 W, (61) 
4tt^ y g - v(q ,q t ) Jq~ 

where g| = lim g0 ^ ±oo g t (g , 0). 

The expression above is valid for single-well potentials in general. Now, let us especialize 
to the potential (|49|). I[q c ] is given by Q55|) , and using ([33]) and (|57D one can write D(q t ,Q) 
as 



D(5 «,e)=(i^r 



A; 2 2fc 2 -l , 
ue H 7^ — E((pe,k) 



k 2 w k 2 

1/2 



cnM 6 dn«e , 2 cnu e snu e ,„* 

+ + 1 - k) — . 62 

sn u@ dn uq 

From ( |53"1) it follows that go — > oo when cn(ue, fc) = 0, which occurs when ue = K(fc), where 



K(k) is the Complete Elliptic Integral of the First Kind. Using (52), this condition can be 
written as an equation in k: 

9 K(k). (63) 



2V2A; 2 - 1 



The graph of f(k) = 2\j2k 2 — 1 K(/c) is plotted in Fig. |2|. It increases monotonically from 
zero (at k = 1/V2) to infinity (as k — > 1), and so for each nonnegative value of G Eq. (|63"D 
has a unique solution, which we denote by k e . Eq. ( p^[ ) then gives the corresponding value 
of g^ (g e = -q£, since U(-q) = U(q)). 



B. Limiting cases of the quadratic approximation 

Expression (EHl) may be used to compute Z 2 (Q) numerically for any value of O. However, 
certain limiting cases may be dealt with analytically, and it is instructive to consider them 
first in order to recover some known results, thus providing a consistency check. These 
limits are: (1) the harmonic oscillator (g — > 0), (2) high temperatures (G — * 0), and (3) low 
temperatures (9 — > oo). 
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1. The harmonic oscillator 



Since V(x) = ^mu 2 x 2 when g = 0, one should obtain the partition function of the 
harmonic oscillator in the limit g — > 0. To see how to arrive at this result starting from the 
expressions derived in this section, we note that in this limit one can perform the integral 
(|61|) using the steepest descent method. For this, we only need the first nontrivial term in 
the series expansion of I[q c ] and D(q t , 0) around q t = 0. 

When \q t \ <C 1, one has k ~ 1 — \ q 2 and k' 2 = \ — k 2 k> ~ q 2 <C 1. Now, we use the 
smallness of k' to derive an approximation for E(^e, k). The first step is to write it in terms 
of k'\ 



E(tp,k) — [ yl — k 2 sin 2 xdx = [ cos x Jl + k' 2 tan 2 x dx. 
Jo Jo 



(64) 



Expanding the integrand to first order in k' , integrating the result, and using the fact that 
siny90 = snwe and cos^e = cnae, we finally obtain 



E(p e ,k) k siiHe + - q 2 



In 



1 + sn uq 
cn-ue 



sn« e 



(65) 



Inserting this result in (|55f) and expanding the square-roots in that equation in powers of 
q t , one obtains 



I[Qc 



Qt 



1 2 1 

- «e In 

2 3 3 



1 + Sn Uq 

cniie 



+ sn uq nc Me 



+ 0(qf). 



(66) 



Now, we use the fact that uq — ► 0/2 and k — > 1 when q t 
cn(u, 1) = sechw, thus finally arriving at 



0, and sn(u, 1) = tanhu and 



I[Qc 



q 2 sinh Q + 0(q 4 t 



(67) 

in (|62l). Then, the first and 



In the case of D(q t , 0) it is enough to take the limit q t 
the last term in square brackets vanish, the second approaches tanh(0/2) [use (p5|)], and 
the third approaches sech 2 (6/2) coth(0/2). Combining these results, one obtains 

D{q u 0) = [Ang tanh(0/2)]- 1 / 2 + 0(q 2 ). (68) 

Inserting (|67|) and (|68 ) in (|6TD yields the desired small-^ limit of the partition function: 

roc p ~{l/2g)q 2 t sinhe 



Z 2 (0) 



Airg tanh(0/2) 



dq t 



2 sinh(0/2)' 



(69) 
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2. High temperatures 

At high temperatures, 6^0 and ( pBD is solved for k& — > l/y/2, and so q@ — ► oo. To 
lowest order in 6, one has cnw e ~ nc«e ~ dnu e « 1, snue ~ «e an d E(<pq, k) « <y? « u e , 
and so the term in curly brackets in fl55j) vanishes in order B: 

/[g c ] = ef/(%) + 0(e 2 ). (ro) 

In -D(g , t, 0), Eq. (0), the third term in square brackets behaves as Uq , while all the others 
behave as Me when — > 0. One has, therefore, 



D(g,,0) = (2^e)- 1 / 2 + O(0 3 / 2 ). (71) 



It follows that 



or, equivalently, 



with V(x) and £7(g) defined in (^Sj) and (£9). This is, clearly, the "classical" limit for the 
partition function with a pre-factor that incorporates quantum fluctuations. 



3. Low temperatures 



At low temperatures, — ► oo and (|63| ) is solved for fee — ► 1. Using the asymptotic 
expansion [Ref. ||40|| , formula 8.113.1] 



1 



K(k) = In (4//c') + - [In (4/Jfe ; ) - 1] fc + 



(74) 



valid for k' = \/i — fc 2 — > 0, it follows from (|63|) that fcg 4e 0//2 in leading order. 
Therefore, 



2/c'| 



4V2e- e / 2 . 



(75) 



Since < g@ -C 1, one has k ~ 1, «e ~ 0/2 and E(<pe, fc) ~ sn(u e , fc) ~ tanh(0/2) 



1 in the whole range of integration. One can also neglect U(q t ) and q^UQ in 



since 



both terms behave as Oe . However, one must be more careful with the combination 
g 2 nc 2 Me (= g 2 ), since it grows very rapidly with q t (in fact diverging when \q t \ — > g@), and 
so cannot be treated formally small quantity. 



Leaving this combination "untouched" in (|55), but making use of the approximations 
listed above, we obtain 



13 



I[q c 



1 2 2 

1 + g % nc Me J 



The analysis of D(q t ,Q) is much simpler. One can simply take the limit 
which was derived under the assumption that g ( < 1. 
Putting all pieces together we finally obtain 



z 2 (e) 



e^oo fit dq 



4 V^9 



t e -i[qc]/g 



(76) 
oo of (|68|), 

(77) 



with g@ and I[q c ] given by (|75|) and (ff6|) , respectively. 



C. Applications 



We shall now apply the quadratic semiclassical approximation to obtain the ground- 
state energy and the curve for the specific heat as a function of temperature. These two 
applications will teach us about the usefulness of the approximation. 

In order to compare (|77|) with the expected low-temperature limit of the partition func- 



tion, Z(Q) 



~ e 



-eeo(s) 



(where Eo(g) = E (g)/fiu; is the dimensionless ground state energy), 



it is convenient to rewrite it in a form in which the 0-dependence can be analyzed more 
easily. This can be done by changing the integration variable back to go- Since q t nctie = qo 
and g@ is the value of q t corresponding to go — ¥ oo, one has 



Z 2 (Q) e ~°° 




exp 




1 2 X3/2 

1+2*° 



1 



(78) 



When ^> 1 it is possible to write an approximate expression for gt(go,0), thus allowing 
to write the integrand in ( [78] ) solely in terms of go and 0. The final result is (see Appendix 
for details) 



Z 2 (0) e ~°° — 



-e/2 



exp < 


r 4 




3/2 

-1 


} 









(79) 



This gives Eo(g) = 1/2, indicating that the quadratic approximation is insufficient to yield 
corrections to the ground state energy of the harmonic oscillator. On the other hand, if one 
recalls that the partition function can be written as 



Z(Q) 



p(0;g,g) dq, 



where 



p(e;g,g) = £e- e£ " 



2 e-oo e _ 6eo 



(80) 



(81) 



is the diagonal element of the density matrix, one may take the square root of the integrand 
in ( [75| ) as an approximation to the (unnormalized) wave function of the ground state. To 
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test the accuracy of this approximation, we have evaluated the expectation values of the 
energy for some values of g and compared them with high precision results found in the 
literature. As Table | shows, the ground state energy computed with this "semiclassical" 
wave function differs from the exact one by less than 1% even for g as large as 2. 

Another concrete problem that can be treated is the calculation of the specific heat of 
the quantum anharmonic oscillator. It can be written in terms of Z(Q) as 



c = e 2 



1 d 2 Z ( 1 dZX 

~zm I ~\zdQ) 



(82) 



This expression was computed using MAPLE for a few values of G and the coupling constant 
value g = 0.3. The result is depicted in Fig. ||, which also exhibits the curve of specific heat 
of the classical anharmonic oscillator (solid line). As expected, the results agree when the 
temperature is sufficiently high, but, in contrast to the classical result, the semiclassical 
approximation is qualitatively correct at low temperatures too, dropping to zero as T — > 0. 

This result, together with the estimate for the ground-state obtained previously, shows 
that the quadratic approximation works very well, being quite accurate at high temperatures, 
and still reliable at lower temperatures. In the next subsection, we will comment on why 
this is so. 



D. Beyond quadratic 



In this subsection, we shall compute the first correction G\ to the quadratic approxima- 
tion. For the sake of comparison, we shall use both the quantum-mechanical (section U. A| ) 
and the field theory approaches (section |II B|) . 

In the quantum-mechanical approach, we want to compute 



Gi(? ,e) 



e 



dd / dr ] G c (0,0;6,r ] )6U[9,r ] }G c (e,r ] ;e,0), 



(83) 



with 8U given by (|58D , and G c constructed from (|56| ) and ([57]), using ( p9[) and (|30"D. Since 
G c is gaussian in 77, and 5U is a polynomial in rj, the integral in rj can be readily performed, 
yielding 



G 1 (q ,Q) = -g 2 



4ngq 2 c (0)Q(0). 



4Q 2 (0) 



(84) 



where the bracket in front can be identified with G c (0, 0; G, 0). 

In the field-theoretic approach, we arrive at the same expression. In fact, the first 
correction corresponds to the m = 1 term in ( |T7| ) which, using fl5Tf), yields 



Gifa,,e) 



e 



de 



(85) 



Using 



we obtain 



15 



Gi(g ,e) = -2 2 G c (o,o ; e,o) f e deg 2 (9,9). (86) 

4 Jo 

((t7 3 ) vanishes, and the factor 3 comes from the three possible pairings of the four 77's in 
(r/ 4 ).) Using (g7|) and (|29|) we reobtain 

Inserting (|86|) in (^l|) and changing the integration variable from qo to (ft gives 

Z(9) = f e dg* D(g t , 9) e- 1 ^ [1 - ^fo, 9) + . . .] , (87) 

J 



where 



ai(g t ,6) = - / d9G 2 (9 } 9). (88) 







Because of the complicated form of Q(9,9), it is not a simple task to compute ai(q t , O). 
However, we can estimate the magnitude of this term without much effort. Indeed, as 
shown in Appendix [|| Q{9, 9) obeys the following inequality: 

0(9,9) < g(9 e 9) (0<9<&). (89) 

Therefore, 

ai(ft,e)<^. (90) 
This shows that (in the case of the quartic anharmonic oscillator) the quadratic approxima- 



tion to the partition function, Eq. ([59]) or (pl|), can be used with confidence whenever the 
condition g9 3 /40 <C 1 is satisfied; this accounts for the numerical agreements obtained in 
the applications of the quadratic approximation. 

A last comment is in order: the next term in the expansion for Z(Q) has a piece with 
a factor g and one with a factor g 2 . The former comes from the product of (1] 6 ) ~ g 3 with 
the overall g~ 2 , whereas the latter involves (r] 8 ) ~ g A . This is another indication that we are 
not dealing with a perturbative series. 



IV. CONCLUSIONS 

In order to understand the nature of the semiclassical series for the partition function, 
it is instructive to look at a diagrammatic expansion of p(Q;q,q), which appears in (|S0D. 
Comparing with (pf), and defining p c (9i, rji, 62, 7/2) = e" 1 ^ 1 ' 92 '^^ 9 G c (9i, 771; 62, 772), we have 



the expansion in Fig. [|, where the dot represents the point q, the squares stand for the 
insertions of [(— l)g -1 5£/], the dashed lines for p c , while the full line represents p(Q,q,q). 
From (^Up and ([H]), we may expand p c itself in terms of the density matrix elements for the 
harmonic oscillator (add and subtract 1 to U"[q c ] in fl2lD ), to obtain Fig. |, where the circles 
represent insertions of {(— l)p -1 [i7"[(7 c ] — 1]}, and the dotted lines stand for the density 
matrix elements of the harmonic oscillator. Therefore, even the first term in our series 
already corresponds to an infinite sum of perturbative diagrams. Analogously, one could use 
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the Feynman diagrams of the field-theoretic description to arrive at the same conclusion. 
This is a clear indication of the nonperturbative nature of our treatment. 

The results of section [TI| can be generalized to higher- dimensional Quantum Statistical 
Mechanics, just as in Quantum Mechanics, where this was accomplished in f26|j27|| . The 
generalization to potentials which allow for more than one classical solution, such as the 
double-well quartic anharmonic oscillator, requires a subtle matching of the series around 
each appropriate saddle-point (i.e., the minima). This is presently under investigation p*[ . 

An extension of our results to field theories is hampered by the fact that we do not 
know how to construct a semiclassical propagator in general. The technical simplifications 
which appear in one dimension cease to exist. However, our methods may still be of use in 
problems where classical solutions have a lot of symmetry (e.g., spherical symmetry) so that 
we can reduce them to effective one-dimensional problems. There are many such examples 
in Physics: instantons, monopoles, vortices and solitons are a few of the backgrounds that 
fall into that category. It is our intention to pursue this line of investigation. 

Finally, we should remark that the field-theoretic treatment of subsection [II B| can be 
used to compute any correlation function of interest. Therefore, a semiclassical series can 
be written down for any physical quantity once it is expressed in terms of the relevant 
correlations. 
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APPENDIX A: 



Here we show that f] b (9), as defined in subsection |II A| , is continuous at 6 = 6/2, and 
that it is also possible to make f)b{0) continuous there by a suitable choice of Q(0). 

To prove that rjb(0) is continuous at 0/2, we have to show that the limits lim^e/2 Vb(@) 
and \im 9 1®/ 2 rjb{9) exist and are equal. However, since rj b (9) is even with respect to reflection 
around 0/2, it suffices to prove the existence of one of them; the existence of the other and 
their equality will then be automatically satisfied. Thus, let us consider the former limit: 



lim T)t,(9) 
ete/2 



lim 

ere/2 



Q(O)g c (0)+g c (0) 



e dff 



o qt 



+ lim v(q, q t ) / 



dq> 



<?o v 3 (q',q t ) 



(Al) 



To compute the above limit, we divide the interval of integration in two subintervals, [go, q] 
and [q,q], with q close to q t but fixed. Then v(q,q t ) times the first integral vanishes when 
q — > q t and, by expanding U(q) around q t , we can approximate v(q,q t ) times the second 
integral by 
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yj2U'(q t )(q-q t ) 



dq' 



h [2U'(q t )(q'-q t )]^- 
Computing the integral and taking the limit, we finally obtain 

1 



lim rib(6) = - 
ete/2 lbK ' U'(q t ) 



(A2) 



(A3) 



Let us now compute the limit 



lim f]h(9) = lim 
eie/2 6»te/2 



6 d6' I 
Q(0) q c (0) + UO) I t^t + 



/o q 2 c (0') UO). 

Using (|T0| ) and ([LTD , and an integration by parts that cancels divergences, yields 

1 /"/•• U"(q) dq 



(A4) 



lim r)b(e) = U'(q t ) 
ete/2 



Q(°) + 7777 



U'(q )v(q ,q t ) J qo [U'{q)f v(q, q t ) J " 



(A5) 



(Although the integrand in ( A~5|) is singular as g — > g t , the singularity is integrable.) Using 
(IE3) we may rewrite ( |A5| ) as 



n 5i# ) = «(o)^,) + ig 



(A6) 



</() 



Since ?)b(#) is odd with respect to reflection around 0/2, lim^e/2 Vb(0) — — h m etG/2 ^fe(^)- 
However, for fjb{d) to be continuous at 0/2, those limits must be equal. This is possible 
only if lim eT e/2^(6 l ) = 0, or 



0(o) 



l 



U'(q t ) [dq t/ 



(A7) 



</0 



APPENDIX B: 



In this appendix, we compute the pre-factor in (|23|). To do so, we remark that ( p0|) is 
the path-integral expression for the Euclidean time evolution operator p c of the quantum- 
mechanical time-dependent problem defined by (note that the role of h is played by g) 



with 



Indeed, we have 



Hm = - 9 ^ 2 + \u"[q c {e)W. 



(Bl) 



(B2) 



G c (9 1 ,ri 1 ;6 2 ,ri2) = {rj 2 \p c {9 2 , 0i)|»7i), 



(B3) 
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so that G c (8i, 771; 62, T]2 ) is a density matrix element. 

Inserting (|23|) into (|B 1| ), and using the fact that I2 satisfies the Hamilton- Jacobi equation 



(B4) 



leads to 



° +1^]g c (9 1 ,0;92,0) = 0, 



(B5) 



\d9 2 2 Q 12 
whose solution is (note that W\ 2 = dVti 2 /d9 2 ) 

G c (9 1 ,0;9 2 ,0) = C(9 1 )n- 2 1/2 . (B6) 
C{8i) can be determined by demanding that we recover the free-particle result as #1 — > 9 2 , 

0i-»02 ro__ i'/i /1 M-l/2 



G c (0i,O;0 2 ,O) D1 ~ 2 [2tT0 (0 2 -0 a ) 
Expanding f2 i2 around 9\, we obtain 



Jll2 — fill + 



90, 



02~- 



$!) + ... = W : 



11 V 172 



0l) + ■ ■ ■ 



(B7) 



(B8) 



Inserting this result in ( |B6| ) and comparing with ( |B7| ) we finally obtain 



G c (9 1 ,0;9 2 ,0) 



W 
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2ng n 12 



1/2 



(B9) 



APPENDIX C: 



Here we provide a simple argument to show why it is usually simpler to write qo{qt, 6) 
rather than q t (qo, 0). In fact, using C = lilt as integration variable in Eq. (0), one has 



6 



loin 2\qt\d( 
1 fimCqt) - U(q t )] 



F qt (q /q t )-F gt (l), 



(CI) 



where F qt (() denotes the primitive of 2\q t \/ y2[U(Cq t ) — U(q t )]. One can then solve ( |CTD for 
q , thus obtaining 



q = q t F- t 1 [e + F qt (l)}. 



(C2) 
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APPENDIX D: 



In this appendix, we derive an approximate expression for qt(qo, 0) valid for ^> 1. We 



start by writing (|T2"D with U (q) given by (f49| ) : 



e = 2V2 ( [ a + / M ) dq = h + J a , (Dl) 

Vl*l / A /(Q 2 -g 2 Ha 2 + ft 2 + 2) 



(q 2 -q?)(q 2 + q? + 2) 

where a is a positive number such that 1 ^> a ^> \q t \ (that both inequalities can be simulta- 
neously satisfied is guaranteed by the fact that \q t \ < g@ ps 4v / 2e _6l//2 when 0>1). Now, 
since both a and |g t | are much smaller than 1, one has q 2 + q 2 + 2 2 for g in the interval 
[|gt|, a], so that 

h - 2 r , dq = 2 cosh- ( ") - 2 In feV (D2) 



9t| yg 2 — g 2 \ l<?< I / \ |9t| . 

where the last (approximate) equality follows from a/ \q t \ being much greater than 1. 

In order to evaluate I 2 , we use the fact that both a and |g | are much greater than \q t \ 
to write 



h „ /« * = 2 In ( M 1 + + « 2 In 2| *' /a . (D3) 

When adding the above results for l\ and I 2 , the dependence on a cancels out (as it should), 
and we are left with the following result: 

0^21*1-^2=1. (D4) 
This can be solved for q t , yielding 



4 on e" e/2 

ft(gb,e)« g ° e . (D5) 



2 

Using this expression to calculate the Jacobian in ([781) one finally obtains (|79|). 

APPENDIX E: 

Let us define the function Qo(9, 6') as 

um = \*j^^ 9 > e , (ei) 

This function satisfies an equation similar to Eq. 
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- ^ Go(0, 0') = 5(9 - 9'), O (O, 9') = o (e, 0') = 0. (E2) 

If we multiply ([10) by Qo(9,9') and (E2) by Q(9,9'), subtract one from the other, and 
integrate the result from 9 = to 9 = O, we obtain 

g (9',9')-g(9',9') = f e d9U"[q c (9)]Oo(0,9 , )g(9,9'). (E3) 

JO 

In the case of the quartic anharmonic oscillator, U"[q c (9)] = 1 + 3q^(9) > and (conse- 
quently) Q(9,9') > for < 9, 9' < Q. Since Qo(9,9') is also nonnegative in this interval, 
Eq. (IEJ) leads to the inequality (BjJ). 
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TABLES 



TABLE I. Ground state energies for different values of g (h = m = u> = 1). 



g E^semiclassical) 3, E (exact) h error(%) 

OA 0.559258 0.559146 0.02 

1.2 0.639765 0.637992 0.28 

2.0 0.701429 0.696176 0.75 

4.0 0.823078 0.803771 2.40 

8.0 1.011928 0.951568 6.34 



a (4>o\H\4>o)/((j)o\(j)o}, where 0o((Zo) is the square root of the integrand in Eq. (f79|), 
b Values quoted from Ref. [Q. 
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FIGURES 
FIG. 1. U{q). 

FIG. 2. Graph of f(k). 

FIG. 3. Specific heat vs. temperature (T = 1/0) for the quantum (diamonds) and classical 
(solid line) anharmonic oscillator, g = 0.3. 

FIG. 4. Diagrammatic expansion for p(Q;q,q). 

FIG. 5. Diagrammatic expansion for p c (Q; qi, q%). 
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